Alternative method for wave propagation analysis within bounded linear media: 

conceptual and practical implications 



O 
O 

(N 



(N 
Oh- 

i 

bp- 

o 

(n : 
> ■ 

(N 



O 

O 



13 



Alberto Lencina/'Q Beatriz Ruiz, 2 and Pablo Vaveliuk 3 

1 Departamento de Fisica, Centro de Ciencias Exatas e da Natureza, 
Universidade Federal da Paraiba, Caixa Postal 5008 CEP 58051-970, Joao Pessoa, Brazil. 
2 Centro de Investigaciones Opticas, cc 124 (1900), La Plata, Argentina. 
s Departamento de Fisica, Universidade Estadual de Feira de Santana, 
Campus Universitdrio BR 166 KM 03, 44031-460 Feira de Santana, Bahia, Brazil. 

(Dated: February 2, 2008) 

This paper uses an alternative approach to study the monochromatic plane wave propagation 
within dielectric and conductor linear media of plane-parallel-faces. This approach introduces the 
time-averaged Poynting vector modulus as field variable. The conceptual implications of this formal- 
ism are that the nonequivalence between the time-averaged Poynting vector and the squared-field 
amplitude modulus is naturally manifested as a consequence of interface effects. Also, two practical 
implications are considered: first, the exact transmittance is compared with that given by the Beer's 
Law, employed commonly in experiments. The departure among them can be significative for certain 
material parameter values. Second, when the exact reflectance is studied for negative permittivity 
slabs, it is show that the high reflectance can be diminished if a small amount of absorption is 
present. 

PACS numbers: 41.20.Jb, 42.25.Bs, 41.85. Ew 



I. INTRODUCTION 

The counter-propagating wave approach is commonly 
used to study the electromagnetic response of spatially 
nondispersive, homogeneous and isotropic plane-parallel- 
faces linear media (Fabry-Perot framework). It consider 
forward and backward monochromatic plane waves to de- 
rive the optical properties as the transmittance and re- 
flectance This method uses the conventional vari- 
ables of electromagnetic fields: moduli and phases of 
these waves which should be found from the Helmholtz 
equation with the corresponding boundary conditions. 
It is well-known that for a single harmonic plane wave 
propagating through an unbounded linear medium, the 
timc-avcragcd Poynting vector modulus is equivalent to 
the squared-field amplitude modulus [2|. This statement 
is generally accepted even if the wave propagation takes 
place in bounded linear media. The counter-propagating 
approach does not point out the possible nonequivalence 
between the above magnitudes since utilizes the ampli- 
tudes and phases as field variables. However, the phase of 
the wave can be replaced by the time-averaged Poynting 
vector modulus as field variable. Then the time-averaged 
Poynting vector and the squared-field amplitude modu- 
lus can be monitored simultaneously within the medium 
so that the conditions that leads to the nonequivalence 
among them will manifest naturally in this equivalent 
frame. 

The so-called S-formalism use the time-averaged Point- 
ing vectos as a field variable. It was recently introduced 
and applied to study the optical response of nonlin- 



mg 

m 



* Corresponding author. Electronic address: agl@fisica.ufpb.br 



ear slabs. The nonequivalence between the time-averaged 
Poynting vector and the squared-field amplitude modu- 
lus was the key to define a nonlinear medium whose non- 
linearity is proportional to the time-averaged Poynting 
vector modulus. Its transmittance was calculated and 
found to differ with that obtained for the Kerr medium, 
whose nonlinearity is proportional to the squared-field 
amplitude modulus. However, the method was not yet 
used to analyze the linear case. Therefore, the aim of 
this paper is to apply the S-Formalism to study the op- 
tical properties of homogeneous, isotropic and spatially 
non-dispersive dielectric as well as Ohmic conductor me- 
dia, inside the Fabry-Perot frame complementing to the 
well-known results of this problem within the counter- 
propagating wave approach. The S-formalism shows, in 
direct form, how this non-equivalence is related to the 
superposition dynamic of opposite traveling waves which 
is only possible in bounded media. Also, the exact trans- 
mittance is compared with that given by the approxi- 
mated Beer's Law, commonly used in experiments, and 
appreciable differences were found for certain material 
parameters values. Moreover, the reflectance of positive 
and negative permittivity media are considered and inter- 
esting results appears for small absorption when negative 
permittivity accounts. 

It must be stressed that the S-formalism does not 
substitute the counter-propagating wave approach be- 
cause the idea of forward and backward waves within the 
medium has a deep physical content. However, the al- 
ternative method adds conceptual and practical implica- 
tions that remains hidden when the conventional scheme 
is utilized. We believe that both methods should be com- 
plementary, helping for a complete physical understand- 
ing on wave propagation in bounded linear media. 

Below, this paper reviews the principal features of the 
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FIG. 1: A harmonic plane wave of amplitude Eq, wave vector 
ko and frequency ui strikes a linear, non-magnetic, homoge- 
neous, isotropic, and spatially non-dispersive plane-parallel- 
faces, to be reflected and transmitted. The field in Region 
II can be represented by a general complex field of real am- 
plitude £ and phase <j>. The Regions I and III constitute, 
for simplicity, the same linear dielectric medium (e.g., air). 
r and t are the reflection and transmission complex coeffi- 
cients, respectively; z — koz is the dimensionless propagation 
coordinate. 



S-Formalism introduced in Ref. |3J as function of the new 
set of field variables. Then, it is applied to derive the 
spatial evolution of both variables for absorbent and non- 
absorbent linear media showing the cases in which the 
equivalence does not hold true, and their causes and con- 
sequences. Finally, practical implications of this features 
are explored: the comparison between the exact trans- 
mittance and the Beer's law, and the analysis of positive 
and negative permittivity media. 



II. S-FORMALISM APPROACH 

When a linear, non-magnetic, homogeneous, isotropic, 
and spatially non-dispersive plane-parallel-faces medium 
of dimensionless thickness d (= kod) is excited perpen- 
dicularly by a plane wave (see Fig. QJ, the S-Formalism 
states that the transmittance and reflectance are re- 
spectively 



R = |rf = 1-5(0), 
T=\t\ 2 = S(d), 



(la) 
(lb) 



where S is the field variable, representing the dimension- 
less time- averaged Poynting vector modulus, given by 



5( - )= <S>_e, = £ 2 # 
Jo dz 



(2) 



where Iq = eoui/(2ko) Eq is the incident intensity with 
eo, the vacuum permittivity, e z the unit vector in the z 
direction, and £ and <j> are the electric field amplitude 
modulus and phase, respectively. S(z) and £(z) evolve 



according to the coupled system Q 



d 2 f 

£ 3 ^+e r £ i -S 2 = Q., 
dz z 

^+a r £ 2 = 0, 
dz 



with boundary conditions 
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J 2=0 
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dz 
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(4b) 
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being e r = e/eo, the relative permittivity, and ay = 
aj (eow), the relative conductivity. From Eq. J2J) it is 
clear that S and £ are equivalents only when cj> is a linear 
function on z. This point is not, in general, emphasized 
in the literature. Thereby, it is mandatory to ask: is 
the equivalence of both magnitudes generalized to cases 
where it could be not longer true?. An example: in Ref. 
0, the energy flux of a TEMoo Gaussian beam propa- 
gating in free-space is calculated from the squared-field 
amplitude modulus instead of the time-averaged Poynt- 
ing vector which could lead to non-physical results. For 
details see Appendix lAl 

Equation i|3b() represents the time-averaged Poynting 
Theorem being ay the responsible for the energy loss in 
the medium. Here we used it explicitly to analyze and 
solve wave propagation problems in bounded media. On 
the other hand, note that, Eqs. (0} are independent of the 
medium properties. Also, the second interface produces 
different values of £ 2 and S at the first interface through 
the spatial dependence of £, as can be seen in Eq. H4afl . 
However at the second interface £ 2 and S are always 
equals with the advertisement that this should not lead 
to false view that they are equivalents. 

Finally note that, the energy conservation is guaran- 
teed thorough the expression: 



\t\ 2 = 1 - 5(0) - S(d) 



(5) 



showing that the reflectance and transmittance values are 
limited by the boundary values of S and enhancing their 
importance in wave propagation problems in bounded 
media. 



III. WAVE PROPAGATION ANALYSIS 

A. Linear dielectric 

The linear dielectric is the simplest nonabsorbent 
medium. It is clear from Eq. I|3bfl that the dimensionless 
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FIG. 2: (Color online) Spatial evolution of the dimensionless 
intensity 5 and squared-field amplitude modulus £ 2 for a lin- 
ear dielectric. Clearly, S is a constant and S 2 is an oscillatory 
function. The parameter values are y/e^ = 1.5, kod — 3tv. 



constant within the medium while the squared-field am- 
plitude is an oscillatory function. Note that £ and S 
are only equals at z r = 1 — 2mir/6 with m e N. Ob- 
serve that when d — rn'ir j '-i/ey (m! e N), 5 — 2m' V, then 
£ (0) = S = 1 attaining their maximum values (T = 1 
and i? = 0). In this case the medium behaves as a <ie- 
lay sheet transmitting all the incident intensity. When 
5 = (2m' + l)ir the field amplitude modulus and the 
temporal-averaged Poynting vector modulus attain their 
minimum values given by £(0) = [l — 4F/(e r — l)]/(l + F) 
and S = 1/(1 + F). Note that, when e r increases, F also 
increases and the minimum of S diminishes. 

On the other hand, when are both magnitudes really 
equivalents? When the squared-field amplitude modulus 
is a constant. In this case Eqs. © relate S and £ by 

S = ^T r £ 2 . (7) 



intensity S is a constant fixed by the boundary condi- 
tions when ay = 0. In this case, the analytical solutions 
for S and £ 2 are 



S = 



1 



£ 2 (Sr) 



1 + Fsin 2 (5/2)' 
l-(4F/(e r -l))sin 2 [(1 



z r )S/2] 



1 + Fsin 2 {5/2) 



(6a) 



(6b) 



where F — (e r — l) 2 /(4e r ) is known as finesse, S = 2dyJlZ 
and z r = z/d. 

Equation l|6a|l is the transmittance result for the lin- 
ear dielectric with Fabry-Perot geometry at normal in- 
cidence, i. e. the well-known Airy-formula p|. The 
solutions © explicitly show that the intensity is non- 
equivalent to the squared-field amplitude modulus. Fig- 
ure |2 points out the nonequi valence. The intensity is a 



This case corresponds to the single plane wave propa- 
gating in an dielectric medium |2|. 

B. Linear absorber 



The linear absorber or Ohmic conductor medium is 
characterized by ay ^ 0. Reference for example, gives 
a rigorous analysis on the wave propagation at Fabry- 
Perot geometry within the counter-propagating wave ap- 
proach. In this frame, additional steps are necessary since 
the time-averaged Poynting vector modulus is not derived 
from the amplitude modulus directly but from the Poynt- 
ing vector definition. In turn, because the S-Formalism 
uses the radiation intensity as a field variable, its values 
at each point within the medium can be directly known 
and compared with the squared-field amplitude modulus 
once the evolution equations are solved. The solutions 
S(z) and £ 2 (z) for the Ohmic conductor are derived from 
S-formalism evolution equations giving (For resolution 
details see Appendix lB|l 



J 



S(z r 



= 2 



a 2 cosh[a_ d(z r — 1)] — a, (£ + 1) sinh[a_ d(z r — 1)] + o?_ cos[a + d(z r 



a_ (£ — 1) sin[« + d(z r — 1)] 



["+ + (£ + !) 2 ] cosh(a_ d) + 2a + (£ + 1) sinh(a_ d) + [a 2 _ 



(£ - l) 2 ] cos (ad) + 2a_ (£ - 1) sin(a d) 



and 



2 (1+0 cosh[a_ d(z r — 1)] — a , sinh[a_ d(z 
£ (Zr) = 4 — : 



1)] + (£ - 1) cos[a + d(5 r - 1)] + a_ sm[a + d(z r - 1)] 



a 2 + (£ + l) 2 ] cosh(a_ d) + 2a + (£ + 1) sinh(a_ d) + [a 2 - (£ - l) 2 ] cos(a + d) + 2a_ (£ - 1) sin(a + d) ' 



(8b) 



where £ = \n c \ and a± = 2 Re 
refraction index given by n 
Eq. © reduces to Eq. © for oy 



2/m{ n c} being n c the complex 
r . It verifies that 



r 



o. 



The solutions iJHJ) are depicted in Fig. |3]fbr different rel- 
ative conductivity values showing that the squared-field 
amplitude modulus could present a markedly oscillatory 
behaviour contrary to the intensity. Figure E^a) shows 
the spatial evolution of both magnitudes for relatively 
small values of ay . Both decrease when the relative spa- 



tial coordinate increases due to the medium absorption. 
The field attenuation is slight for small conductivity val- 
ues which implies an intense contribution to the total 
field of the wave returning from the second interface pro- 
ducing a strong wave superposition. On the other hand, 
Fig. |2Ib) was depicted for intermediary ay value. The 
field attenuation is sufficiently strong and the field ampli- 
tude modulus acquires a quasi-negligible value at inter- 
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face z r = 1. However, this still produces a wave superpo- 
sition so that £ 2 shows a slightly oscillatory decreasing 
behaviour. Thereby, the intensity and the squared-field 
amplitude modulus are yet nonequivalents. Finally, Fig. 
|3fc) was depicted for a relatively high value of ay. In this 
case, the field amplitude modulus is quite attenuated be- 
fore to reach the interface z r = 1. Then, the medium 
behaves as unbounded and both, the intensity and the 
squared-field amplitude, become completely equivalents 
verifying Q 

S{z r ) = ^±£ i {z r ). (9) 

In summary, Fig. [3]shows that the strength of wave su- 
perposition dynamic is directly related with the effect of 
the back interface which turns out in the nonequivalence 
of S and £ 2 . 



IV. BEER'S LAW COMPARISON 

The Beer's Law is commonly used to calculate the ab- 
sorption coefficient by measuring the transmittance when 
the optical absorption within the medium is accounted. 
This law gives the intensity of the wave at interface z r = 1 
when the second interface effects arc neglected 6] . There- 
fore, this Law is an approximate result since it consider 
only the first interface (only one boundary condition) to 
calculate the energy flux. The dimensionless intensity 
given by Beer's Law Ib attenuates exponentially as 

I B (z r ) = S(0)e- Q - 5 ", (10) 

where a_ characterize the absorptive medium properties 
(note that was used S(0) instead of Iq supposing that 
the intensity at z r — is know and it lead to dimin- 
ish the differences with the exact result). Equation i|10|) 
was also depicted in Fig. [21 with the aim to compare 
the approximate Beer's Law intensity and the true in- 
tensity within the medium. It happens that Ib{z t ) can 
departure from S(z r ) for a wide range of usual parame- 
ters. Thereby, the exact transmittance S(l) could differ 
from the approximated transmittance value used 
often in experiments. It is clear that a difference appears 
when the second interface mediates in the wave dynamic, 
as Fig. Ola) shows. When the effect of that interface is 
either quasi-negligible [Fig. ISJb)] or the medium can be 
considered as unbounded [Fig. Etc)], both Ib(zv) and 
S(z r ) do not present differences. Thereby, the validity 
of this approximation must be carefully tested for each 
particular problem since the difference S(l) — -fs(l) de- 
pends on material parameters: e r , ay and d. In Fig. 0] 
the percentaged difference — 5"(1) is shown. As 

it can see, there exist cases where the difference of both 
transmittances is approximately 10%. This major depar- 
ture appears when the back interface effect is relevant 
i. e., when the wave superposition dynamic plays an 
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FIG. 3: (Color online) Spatial evolution of the dimensionless 
intensity S, squared-field amplitude modulus £ 2 , and Beer's 
Law intensity Ib for a linear absorber with different relative 
conductivity values: (a) oy = 0.1, (b) 1, (c) 10, with = 
2.5, and d = 2tt. Clearly, the non-equivalence between S and 
£ 2 is directly related to the effect of the second interface. 

important role which happens only for certain particu- 
lar values of (e r ,oy). Notice that, the maximum differ- 
ence appears for high permittivity values and low, but 
nonzero, absorption. 

V. REFLECTANCE OF POSITIVE AND 
NEGATIVE PERMITTIVITY MEDIA 

Negative permittivity media have attracted the atten- 
tion of scientist some years ago because was indicated a 
way to build them its |j| and because was open the pos- 
sibility of have others media like the Left-handed Media 
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FIG. 4: (Color online) Percentaged difference Jb(1) - S(l) 
against e r and ay for d — 2n. 

8] leading to a intense research in this area 0, 0, • 
However, always the reflectance of negative permittivity 
media was analyzed in the context of the effective permit- 
tivity of a given microscopic configuration 0, [Hi • Here 
we done a completely macroscopic analysis of the opti- 
cal properties of positive and negative permittivity media 
without enter in the details of the microscopic configura- 
tion. 

In Fig. El (a) the reflectance for both, positive and 
negative values of the relative permittivity is shown for 
several orders of magnitude of ay. By observing this fig- 
ure, we found three regions well differentiated: (1) e r > 
and ay < 1; (2) e r < and ay < 1; (3)oy > 1. Region (1) 
is characterized by an oscillating low reflectance whose 
envelope grows as far as e r increases. This is a well-know 
behavior the reflectance of a low loss Fabry-Perot. 
On the other hand, region (2) has a well particular de- 
pendence on e r and ay: the reflectance is close to unity. 
Finally, region (3) has a uniform reflectance as function 
of e r that grows monotonously as ay increases, such as in 
the limit ay 3> 1, R ~ 1 no matter the e r -value. 

For a deeper understanding of the reflectance behavior 
we analyze the a± dependence on e r and ay. Figures 
03(b) and [SJ (c) depicts a+ and a_, respectively, as func- 
tion of e r and oy. Also, in Table I the limiting values 
for ay <C of a± and R are shown, to helping in 

the analysis. Note that, a± represent the effective per- 
mittivity and conductivity in all the regions. Therefore, 
an analysis of these magnitudes becomes important for 
a better physical understanding on the wave dynamic 
propagation. In fact, these magnitudes determine the 
characteristics of the three region mentioned above, since 
their values changes dramatically in those. In region (1), 
a + 7^ and a_ = 0; in region (2) a + = and a_ 7^ 0; 
and in region (3) a+ ~ a_ =/= 0. Then, the low and 
oscillatory reflection correspond to the region where a+ 
predominates, i.e. the medium effectively behaves as a 
dielectric where the oscillations are produced by the wave 
superposition as result of the existence of the second in- 
terface summarizing in an Airy-type behavior, as Table I 



TABLE I: Limiting values for a± and R for the three regions 
mentioned in the text. 
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shows. On the other hand, it is well-known that the high 
reflectance is commonly associated with high values of 
the conductivity, i.e for good conductors. However, high 
reflectance could occur even for low ay when the permit- 
tivity is negative. In this case, the response is dominated 
by a_. The reflectance is given by a hyperbolic- Airy- 
type function as Table I show. Because F' sinh(5/2) 3> 1 
for almost all the e r -values in this region, then R rj 1 
and the medium posses mirror-like properties. For region 
(3), a± strongly depend on ay-values. The medium re- 
sponse is completely dominated by ay with a reflectance 
that asymptotically reaches the unity with a moderately 
slow rate. Indeed, the asymptotic high reflection region 
is dictated by a unique parameter a ~ a+ w a_. Sum- 
marizing, for region (1) is obtained a periodic function 
of e r for the reflectance, whereas for region (2) it con- 
tains sine hyperbolic functions that monotonously grow 
as e r decreases achieving rapidly the unity. Note that, al- 
though each reflectance was obtained by calculating the 
appropriate limit of Eq. JSJ , the reflectance in region (2) 
could be obtained from that of region (1) by changing 
e r — > — e r what allows us to say that the propagation pa- 
rameter a+ "becomes" imaginary when the permittivity 
takes negative values doing that the medium behaves as 
a high absorbing one. 

To end this section we stand out others interesting 
results: for e r = 1, R = because the medium is 
"matched" with the vacuum and all the light is trans- 
mitted; For e r = and a r = 0, R = (1 + 4/d)~ 1 , showing 
that also in this "quasi-nihility" |l2( some quantity of 
electromagnetic energy can propagate. And finally, note 
that for e r < and ay ~ 1 there exist a region where 
R ~ 0.5 showing that not all the negative permittivity 
media has a high reflectance, but those with oy <C 1. 
Similar results were observer for low loss negative per- 
meability media . 



VI. CONCLUSIONS 

The S-Formalism was used to study the wave propa- 
gation in bounded linear media. It utilizes the set (5, £ ) 
as field variables instead of the conventional set (cf>,£). 
Then, the electromagnetic energy flux can be directly 
known inside the medium. The approach stresses that, in 
general, the time-averaged Poynting vector modulus and 
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FIG. 5: (Color online) Reflectance (a), Logioa+ (b), and LogioQ- (c) against e r and Logiooy. For reflectance d — 2ir. 



the squared- field amplitude modulus are non-equivalents. 
The analysis clearly shows that the non-equivalence of S 
and £ 2 takes place in bounded media being consequence 
of the back interface. The role of the latter is responsible 
for the superposition dynamic between the forward and 
backward waves which causes an oscillatory behaviour 
of the squared-field amplitude modulus contrary to the 
intensity. When the second interface can be neglected, 
the medium can be treated as unbounded with a single 
wave propagating such that S and E 2 are equivalents. 
The analysis shows that the usual Beer's Law approxi- 
mated intensity, employed to calculate the transmittance 
in usual experiments, could departure significatively from 
exact transmittance and this difference is also produced 
by the effect of the second interface. Thereby, the valid- 
ity of this approximation should be rigourously tested for 
each particular problem. Moreover, positive and negative 
permittivity media were analyzed and found different be- 
haviors accordingly to the permittivity and conductivity 
values. Was observed that, for low absorption, negative 
and positive permittivity media have a well differentiated 
reflectivity (being unity for the former and a low oscil- 
lating one for the later) whereas on the contrary case the 
reflectance does not depend on the permittivity values. 
For intermediate (moderate) absorption, the reflectance 
can take values around 0.5, showing that not all negative 
permittivity media have high reflectance. 

Both methods, the conventional and the S-formalism, 
should be used complementarily what could help for a 
complete physical understanding on wave propagation in 
bounded linear media. The results for finite bandwidth 
waves can be easily obtained from the monochromatic 
ones by i nteg rating all the contributions for the trans- 
mittance 13] . 
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APPENDIX A: TRUE POWER OF A TEM 0() 
PARAXIAL GAUSSIAN BEAM 

The free-space propagation of a TEMoo Gaussian beam 
is a good example on two fundamental aspects: first, 
because the misconception about the equivalence of the 
squared-field amplitude modulus and the time-averaged 
Poynting vector exists in the literature and second, be- 
cause by applying this misconception could lead to erro- 
neous physical results. 

The theoretical framework on Gaussian beams within 
the transverse field-paraxial approximation is well-known 
H 0. The structure of a TEMoo held amplitude, 
E(x,y,z), propagating in free-space along_z-axis with 
wave vector modulus k (towards +z) is f 



E 



WQ 

w (z) 



exp 



y 



w 2 (z) 



exp [ikz] 



x exp 



-i tan 



exp 



,,21 



ik 



2R(z) 



(Al) 



here by simplicity we use dimensional coordinates. The 
Gaussian TEMoo beam wavefront is perfectly flat at 
z = acquiring, thus, curvature and begin spreading 
in accordance with 



R{z) = 



z + 



w 2 (z) 



1+ - 



,(A2) 



where z$ = £;u>q/2, and z is the distance propagated from 
the plane z = where the wavefront is flat being wq the 
minimum spot size often called the waist radius. The pa- 
rameters w(z) and R{z) are the spot size and the wave- 
front radius of curvature respectively after the wave has 
propagated a distance z 0, |(|- In Ref. 4, it has been 
emphasized that the energy per unit time crossing an in- 
finite transverse plane (z = cte) should be a constant to 
meet the energy conservation for a losses medium. This 
calculus has been done considering equivalents the time- 
averaged Poynting vector and the squared-field ampli- 
tude modulus since 
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FIG. 6: Both, P m and Pg& Gaussian beam powers against 
waist radius for Eq = 1.2 x 10 s V/m. The inset depicts the 
ratio among both powers as function of dimensionless param- 
eter kwo 



was used giving the total power 0: 



p -^ £ 2 

_ 4fc 



(A3b) 



which effectively is a constant. However, P m does not 
represent the truly total power as it will be shown in the 
following. Necessarily, the power should be calculated 
from the time-averaged Poynting vector as 



Pt 



Gb 



(S z ) dxdy, 



(A4a) 



— oo — oo 



where (S z ) is the z-component of the time-averaged 
Poynting vector given by Eq. @. Calculating (|A4a|) 
it has 



we uJ 2 W 2 



1 



(A4b) 



so that the power Pcb is different of the power P m ac- 
counted in the literature. 

Figure depicts Pcb and P m against the minimum 
spot size simulating an Argon laser beam. In most of the 
experimental cases, the beam wavelength is much lesser 
than that the beam waist radius satisfying kwo ^> 1. In 
these cases Pcb ~ Pm, as the inset of Fig. [fj] shows. It is 
clear that both power curves disjoin only from kwo < 10. 
This condition could be satisfied by ultra-focused beams 
[l4T| . In particular, note that, kw < 1 leads to Pcb < 
(switch in the beam energy flux direction) which, of 
course, is a non-physical result. This unreal value of the 
total power Pcb would be consequence of the fact that 
the transverse field-paraxial Gaussian beam approach are 
not longer valid. Thereby, the power curve Pcb could 
be useful as measure parameter to evaluate when this 
theory can be used: when the Pcfc-value is far of the 
values equal and lesser than zero, i. e. when Pq^ 
P m . For the example placed in Fig. this happens 
from wq ~ 300 nm. On the contrary, it is found that 



always P m > 0, which could lead to a misconception that 
this value represents the real beam power. The power 
curve P m does not indicate when the theory fails. In 
summary, the departure of Pcb from P m could indicate 
the validity degree of the traverse field-paraxial Gaussian 
beam approach. 

This simples example shows that the misconception 
about the equivalence of that magnitudes exists. In the 
literature, often, \E\ 2 is used to calculate the intensity 
which could lead to erroneous physical results. 



APPENDIX B: RESOLUTION OF EQS. © 

Here is show the detailed resolution of Eqs. Thus, 
the dielectric case is easily obtained making oy = 0. 
Equations © can be written as: 



n d 2 u ( du\ 2 2 jo2 

2 d^ u -{dE) + 4 ^ 2 = 4 ^ 

dS 

— = -a r u, 
dz 

whit u = £ 2 and the boundary conditions reads 



du 
dz 



. o 1 du 
(u + S) 2 + -— -4m 
Adz 



= 0, 
0. 
0. 



(Bla) 
(Bib) 

(B2a) 
(B2b) 
(B2c) 



z—d 



By the homogeneity of the Eqs. \B1\) . we write down 
the following ansatz: 



u = A exp (az) , 
v = B exp {az ) . 

Replacing in Eq. HB1|) . it has 

a 2 = idi, 



(B3a) 
(B3b) 

(B4) 



Then, the general solution can be written (by the homo- 
geneity of Eqs. f \Bl\) ) as a linear combination of the four 
possible values of a: 



u = Aiexp (ot-Z) + A 2 exp (— oc-Z) + 

+A 3 exp (ia + z) + A^exp (~ia + z) 



a+ot- A\ A 2 . 
exp{oi-Z)-\ exp(— a_z) - 

+i — -exp (ikoa+z) — i — -exp (— ia+z) 
a + a+ 



(B5a) 



(B5b) 



The boundary conditions fix three of the four constants 
appearing in the solution, the fourth is fixed by auto- 
consistency. Replacing Eq. l|B5a(l in Eq. i|Bla|) , it results 



AiA 2 



Ax A 



3^4- 



(B6) 



8 



Applying the boundary conditions at z = d and defining 
B = A 1 /A 3 , C = A 2 /A 3 , results BC = A 4 /A 3 , has 



B 



C = 



2i£ — (a_ + ia + ) 
2i£+ (a_ +ia + ) 

2i£ — (a_ — ia + ) 
2i£, + (a_ — ia+) 



exp 



exp 



(— a-+ia + )d , (B7a) 



(a_ + ia + ) d , (B7b) 



and with the boundary condition at z = 

A 3 =4[B(l + £-a + ) + C(l + £ + "+) + 

+ J BC(l-^-ia_) + (l-e + ja-)] _1 . (B8) 



I 



By replacing this values in Eq. IjBSafl . it gives 



u = 4 



(1 + £ — a + ) exp [a_(z r — 1)] + (1 + £ + a + ) exp [— a_ (S r - 1)] + 
(1 + £- a ) 2 exp[-a_d] + (! + £ + a,) 2 exp[a_J]- 



-(^ — 1 + ia_ ) exp [— ia + (z r — 1)1 + — 1 — ia_ ) exp [za + (z r — 1)] 
— (£ — 1 + ia_) 2 exp[ia + (i] — (£ — 1 — ia_) 2 exp[— ia + d] 



(B9a) 



and 



u = 2 



(1 + £ + exp [— a_ (z r — 1)] — (1 + £ — a + )a + exp [a_ (z r — 1) 

(1 + £- a ) 2 exp[-a„J] + (1 + £ + a + ) 2 exp[a_d\- 



— 1 — )ck_ exp [za + (5 r — 1)] — — 1 + ia_ )a_ exp [— za + (z r — 1)] 
— (£ — 1 + ia_ ) 2 exp[ia + ii] — (£ — 1 — ia_ ) 2 exp[— ia,a| 



(B9b) 



From this equations, Eqs. (JHJ can be obtained rearranging the exponential terms. 
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